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^ ; ABSTRACT 

rr . We report millimeter-wavelength observations of the Sunyaev-Zel'dovich (S-Z) effect 

^_. I in two distant galaxy clusters. A relativistically correct analysis of the S-Z data is com- 

^ ' bined with the results of X-ray observations to determine the radial peculiar velocities 

(vr) of the clusters. We observed Abell 2163 {z = .201) in three mm-wavelength bands 
^ . centered at 2.1, 1.4, and 1.1mm. We report a significant detection of the thermal 

cn I component of the S-Z effect seen as both a decrement in the brightness of the CMB at 

^ ■ 2.1mm, and as an increment at 1.1mm. Including uncertainties due to the calibration 

^ ■ of the instrument, distribution and temperature of the IC gas, and astrophysical confu- 

O ! sion, a simultaneous fit to the data in all three bands gives Vr = +4:90llll™ kms~^ at 68% 

confidence. We observed Abell 1689 {z = .181) in the 2.1 and 1.4 mm bands. Including 
the same detailed accounting of uncertainty, a simultaneous fit to the data in both bands 






% 



Q_^^ gives Vr = +170_63okms . The limits on the peculiar velocities of A2163 and A1689 

X ', correspond to deviations from the uniform Hubble flow of < 2 — 3%. 

I : «.. ...... — .: ..™.o. - ... ._ .„ - ..... 

clusters: individual (Abell 2163, Abell 1689) — X-rays: galaxies — large-scale structure 
of the universe 
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1. Introduction 

Inhomogeneities in the large-scale mass distribu- 
tion of the universe produce, through their gravi- 
tational interaction, a corresponding velocity field. 
Measurements of this velocity field can place use- 
ful constraints on cosmological mode ls, includin g 
the mass density of the universe (e.g., Dekel 1994 ). 



Galaxy clusters have been shown to be effective trac- 
ers of the large-scale velocity field of the universe 
( Bahcall et al. 1994 ; Gramann et al. 1995 ). Stan- 



dard methods for determining the radial peculiar ve- 
locities [vr) of clusters difference the velocity deter- 
mined from the redshift [z) with that expected for 
the uniform Hubble flow, 



Vr 



<il 



Qoz + {qo - 1) i^V^ 



z + l-l) -HodL, 

(1) 

where Hq is the Hubble constant and qo is the de- 
celeration parameter. The luminosity distance to the 
cluster (di) is typically determined with an empiri- 
cal relationship relating the brightness of a standard 
candle to its distance. In general, the uncertainties 
in the distance indicators, and therefore, the cluster 
peculiar velocities; increase linearl^r with distance 



The interaction of the Cosmic Microwave Cack - 



ground (CMB) with hot intracluster (IC) gas bound 
to clusters of galaxies provides a method to directly 
det ermine the radial peculiar vclocitica of galaxy clua — 



teri: — Compton acattering of the CMB by hot IC 
gas - the Sunyaev-Zel'dovich (S-Z) effect - gives rise 
to an observable distortion of the CMB spectrum. 
For a general review of the S-Z effect, see Rephaeli 
(1995a). The change in CMB intensity caused by 
the scattering has a thermal component due to the 



random motions of the scattering electrons (Sunyaev 



fe Zel'dovich 1972), and a kinematic component due 



to the bulk peculiar velo city of the cluster gas ( Sun- 
^aev fc Zel'dovich 1980| ). The S-Z thermal compo- 
nent produces a decrement in CMB intensity at low 
frequency and an increment at high frequency, while 
the kinematic component appears as an increment 
or decrement at all frequencies. In Figure |l|, we 
show the change in the CMB intensity for the two 
components of the S-Z effect as a function of fre- 
quency. Accurate measurements at two or more mil- 
limeter/submillimeter wavelengths can, in principle, 
separate the two components of the S-Z effect. 

When combined with X-ray measurements, the 
amplitude of the thermal component can be used to 
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Fig. 1. — Change in CMB Intensity due to the two 
components of the S- Z effect. The non-relativistic 
approximation to the thermal component (solid line) 
is plotted for y = 10^"* and the kinematic component 
(dashed line) is plotted for {vr — ±1000 kms~ , r = 
0.01). The average transmissions of the three SuZIE 
spectral bands (normalized to one) are superimposed. 



determine Hq (Cavaliere et al. 1979; Birkinshaw ei 



1991). The ratio of the kinematic and thermal 



component intensities can be combined with the IC 
gas temperature to determine the the radial com- 
ponent of the cluster peculiar velocity (3unyaev & 
Zel'dovich 1980; Rephaeli & Lahav 1991). Because 



the surface brightness of the S-Z effect is indepen- 
dent of redshift and the derived peculiar velocity is 
weakly dependent on IC gas distribution and temper- 
ature, accurate peculiar velocities can be determined 
for distant clusters which, in general, have poor X-ray 
data. 

We have developed an instrument and observing 
strategy which achieves high sensitivity to the S-Z ef- 
fect in three mm- wavelength bands while minimizing 
systematic errors ( Holzapfel et al. 1997a ). This pa- 
per describes the observations and analysis used to 
place upper limits on the peculiar velocities of two 
distant galaxy clusters. In Section ^, we discuss the 
spectra of the S-Z effect and the determination of pe- 
culiar velocities. The instrument used to make the 
S-Z observations and its calibration are described in 
Section p. . Analysis of the S-Z data is described in 
Section ^. The results of the observations of A2163 
and A1689 are described in Sections ^ and |6]. Astro- 
physical confusion of the S-Z effect, including the con- 



tribution of primary CMB anisotropies, is discussed 
in Section Uj. We summarize our results and conclu- 
sions in Section ^. 

2. Theory 

2.1. Relativistic Comptonization in Clusters 

The Sunyaev & Zel'dovich (1972) treatment of 
CMB Comptonization, which is commonly used in the 
analysis of the S-Z effect in clusters, is valid only at 
low gas temperatures. The determination of accurate 
cluster peculiar velocities using the S-Z effect requires 
the use of the relativistically correct treatment of the 
Comptonization. We express the relativistically cor- 
rect intensity difference as 



AJt = /, 



J^ix,Te)dT, 



(2) 



where Iq = 2(kTQf/{hcf, Tq is the CMB tempera- 
ture, dr = UpCTTdl, and the integral is along the line 
of sight through the cluster. The spectrum is given 

by 



*(:J, Te) 



['^{x,T,)~l] 



(3) 



(e^ - 1) 

where $(a; , Te) is a three dimensional integral fully 
specified in Rephaeli (1995b). To facilitate compar- 
ison of our results with low frequency observations, 
we determine an effective Comptonization in terms of 
the measured surface brightness. 
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As we will demonstrate in Section 3.3.2., neglect 



ing relativistic corrections would lead to considerable 
errors in the cluster peculiar velocity. All results are 
calculated using the relativistically correct treatment 
of the Comptonization unless otherwise stated. 

2.2. Peculiar Velocity from S-Z and X-Ray 
Measurements 

If the IC gas has a bulk velocity (vp) with respect 
to the CMB frame, there is an additional kinematic 
component to the S-Z effect. 



AIk = —lo h{x) / UeCTT • dl , 



(5) 



where the spectrum is given by /i (a;) = x'*e^/(e^ — 1)^, 
identical to that for a change in the CMB tempera- 
ture. 



Following Birkinshaw, Hughes & Arnaud (1991) 
(BHA), we express the IC gas temperature and den- 
sity in terms of the product of a reference value and a 
dimensionless form factor (BHA equations 3.3 - 3.4): 



rie = rieofniO, 0, ^), 



(6) 
(7) 



The angle 9 is measured from a reference line of sight, 
(j) is the azimuthal angle about that line of sight, and ^ 
is the angular distance along the line of sight. We ex- 
tend this treatment to include a term describing the 
temperature dependence of relativistic Comptoniza- 
tion, 

*(x, Te) = *o(a^, T,o)h[fT{e, 0, 0]- (8) 

The expressions for the intensity of the S-Z compo- 
nents, equations ^ and |5|, then become: 
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The structural information for the cluster is contained 
in the angular form factors: 
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(11) 



(12) 



The peculiar velocity can be expressed in terms of 
the intensity ratio of the two components of the S-Z 
effect, 



c ~ AlTh{x)eK 



(13) 



In general, measurements are made over finite 
bands in frequency which contain contributions from 
both components of the S-Z effect. The intensities 
of the two components of the S-Z effect can be ex- 
pressed as linear combinations of the intensities in 
the measured bands. In the case where the IC gas is 
isothermal, the angular scale factors Qt and Qk are 
identical and results for Vr depend only on the ratios 
of the measured band intensities. 



Observations of the S-Z Effect 



3.1. Instrument 



The S-Z observations were made using the 
Sunyaev-Zel'dovich Infrared Experiment (SuZIE) 
bolometer array at the Caltech Submillimeter Obser- 
vatory (CSO) on Mauna Kea. A detailed description 
of the SuZIE instrument has been presented elsewhere 
([Holzapfel et al. 1997a|). SuZIE is a 2 x 3 array of 



300 mK bolometric detectors optimized for observa- 
tions of the S-Z effect in distant (z > .1) clusters of 
galaxies. The array consists of two rows separated 
by 2.2' in declination; each row consists of three co- 
linear array elements separated by 2.3' in right ascen- 
sion (RA). Each array element produces a beam on 
the sky which is approximately 1.7' FWHM. 

The spectral responses of the array elements are 
determined by a common set of metal-mesh filters. 
The filters can be changed for observations in three 
mm-wavelength passbands. The 2.1 and 1.1 mm filter 
passbands are designed to maximize the ratio of the 
S-Z thermal component signal to the sum of the at- 
mosphere and detector noise. The 1.4 mm filter band 
is designed to have high sensitivity to the S-Z kine- 
matic component while constraining the net thermal 
component contribution to be near zero. Due to rel- 
ativistic corrections, this cannot be done for all IC 
gas temperatures; in general, the 1.4 mm band con- 
tains some residual thermal component signal (see 
Section 3.3.2.). The average transmission of each of 
the three SuZIE pass-bands are shown in Figure El 

Array elements within each of the two rows are 
electronically differenced by placing pairs of detec- 
tors in AC biased bridge circuits. The output of each 
bridge is synchronously demodulated to produce a 
stable DC voltage proportional to the instantaneous 
brightness difference between the two array elements. 
The bolometer differences strongly reject signals com- 
mon to both elements, such as fluctuations in the tem- 
perature of the 300 mK heat sink and atmospheric 
emission. Each row provides three differences, two 
with beam throws of 2.3' and one with 4.6'. 

3.2. Scan Strategy 

For observations of the S-Z effect, the SuZIE array 
is used in drift scanning mode. The telescope first 
tracks a position leading the source by a right ascen- 
sion offset (RAO). To begin the scan, the telescope 
stops tracking and is fixed in place; the rotation of 



the earth drifts the source across the array of detec- 
tors. The scans are roughly 30' in length and last 
~ 120/ (cos^) sec for low declination (6) clusters. Be- 
tween scans, the array is rotated about the optical 
axis so that the two rows of the array are kept paral- 
lel to the direction of the scan. Keeping the telescope 
fixed while taking data eliminates signals due to mod- 
ulation of the telescope's side-lobes and microphonic 
response of the detectors. After a scan, the telescope 
tracks a new position leading the source by a differ- 
ent RAO and a new scan is begun. This cycle, which 
alternates the RA O be tween scans is repeated indef- 
initely. In Section 4.6. , we use the alternation of the 
RAO between scans to test for an instrumental base- 
line. 

3.3. Calibration 

3.3.1. Flux Calibration 

Observations of planets are used to map the beam 
shapes of the instrument and calibrate the responsiv- 
ity of the detector elements. In April 1993 and 1994, 
Uranus was used to calibrate the 2.1 mm observations. 
We used scans of Uranus to calibrate the 1.4 mm band 
in May 1994. In May 1993, observations of Mars were 
used to calibrate 1.1 mm band observations. Rotation 
of the array about the optical axis results in small re- 
producible changes in the instrument's beam-shapes. 
For each band, we mapped the beam shapes over the 
range of dewar rotation angles for which we observed 
the source; the calibrations computed for the range 
of rotation angles at which we observed the sources 
change by less than 5%. 

The brightness of Uranus is found from a third or- 
der polynomial fit to the measured brightness tem- 
perature as a function of wavelength ([Griffin and Or- 
ton 199tj ). The brightness temperature of Mars is 



taken from Orton et al. (1986). The authors assign 
±6% uncertainty to the brightness of Uranus, most 
of which arises from the ±5% uncertainty in the ab- 
solute brightness of Mars from which Uranus is cali- 
brated. Combining the uncertainties in the measured 
beam-shapes and absolute brightness of the planetary 
calibrator, the uncertainty in the calibration of the in- 
strument to surface brightness in each of the spectral 
bands, is estimated to be ~ ±8%. 

We use these results to calculate the contribution 
of the uncertainty in flux calibration to the uncer- 
tainty in the peculiar velocity determined from mea- 
surements in two spectral bands. Because the pe- 



culiar velocity can be expressed as a ratio of mea- 
sured brightnesses, a calibration error which is com- 
mon across the bands, such as the uncertainty in 
the absolute calibration of the planetary calibrator 
is likely to contribute, makes no contribution to the 
uncertainty. We have calculated the flux calibration 
error assuming the worst (and unlikely) case that the 
< 5% uncertainty in the measured beam-shapes is an- 
ticorrelated between bands. For the combination of 
the 2.1 and 1.1 mm data, we determine the contri- 
bution to the uncertainty in the peculiar velocity to 
be Aw^ ^ ±140(/cTe/10keV)kms"^ For the combi- 
nation of the 2.1 and 1.4 mm data, the contribution 
to the uncertainty is Avr < ±6(fcTe/10keV) kms~ . 
The comparatively large uncertainty determined for 
the 2.1 and 1.1mm data is because, to determine 
Vr, we are (in effect) differencing two large thermal 
component signals, each with a proportional flux cal- 
ibration uncertainty. For the combination of the 2.1 
and 1.4 mm data, the uncertainty is roughly propor- 
tional to the small signal in the 1.4 mm band. For 
a peculiar velocity determination using the 2.1, 1.4, 
and 1.1mm bands, the uncertainty in the flux cal- 
ibration contributes an uncertainty Avr ^ ±(6 — 
140) (Te/lOkeV) kms~^, depending on the relative 
significance of the data in the three bands. In our 
accounting of the uncertainties we assume the worst 
case, and use the upper limit. In this analysis, we 
have assumed that 5% of Vr is a negligible contribu- 
tion to the total uncertainty, true for the results of 
this work. 

3.3.2. Spectral Calibration 

The spectral responses of the array elements, in- 
cluding detailed checks for out of band leaks, are mea- 
sured using a Fourier transform spectrometer. Be- 
cause the beams of the array look through different 
portions of the filter at different angles, there are 
slight differences between their spectral bands. Small 
deviations of the 1.4 mm band from the null frequency 
of the thermal effect produce large offsets in the de- 
rived peculiar velocities. To demonstrate the size of 
this effect, we calculate the peculiar velocity that one 
would determine by assuming all the signal in the 
1.4 mm band is due to the kinematic effect. 



AVr ^ ^dv^{u,Te)f{v) 

c J dv h^h') f {i>) 



(14) 



In the non-relatistic approximation, the thermal ef- 
fect signal in the average 1.4 mm band is nearly zero. 



In Table |l|, the scatter of the non-relativistic results 
about zero are due to the array element differences. 
The results are listed for Tg = lOkeV; the offsets 
scale linearly with the assumed Tg. These significant 
differences in the spectral response of the array ele- 
ments are taken into account by individually calibrat- 
ing each of them to the two components of the S-Z 
effect. 

Relativistic corrections to the S-Z spectrum lead 
to significant thermal component signal in the 1.4 mm 
band. In Table El we show the error in peculiar ve- 
locity that would be incurred by ignoring relativistic 
corrections to the S-Z effect for the three principal 
gas temperatures considered in this work. In the non- 
relativistic approximation, the average 1.4 mm band 
is centered at the null of the S-Z thermal effect. Rela- 
tivistic corrections shift the null to higher frequencies; 
the change in freque ncy is linearly d ependent on the 
IC gas temperature ( Rephaeli 1995a ). The S-Z spec- 
tra, as sampled by the narrow 1.4 mm band (See Fig- 
ure y), increases nearly linearly with frequency. The 
amplitude of the thermal component also increases 
nearly linearly with increasing IC gas temperature. 
Combining these effects, we find that neglecting rel- 
ativistic corrections leads to an error in the peculiar 
velocity determined form the average 1.4 mm band of 
Awr « -h360(/cTe/10keV)2. This is corrected for by 
calibrating with the correct relativistic spectra. 

There is uncertainty in the results of the laboratory 
spectroscopy. From repeated spectroscopy runs, we 
conservatively estimate the uncertainty in the array 
element central frequencies to be < ±1%. When the 
peculiar velocity is determined from the combination 
of the 2.1 and 1.1mm band data, the uncertainty in 
Vr is ±137(fcre/10keV)kms"\ The steeper slope of 
the thermal spectra near the null leads to a slightly 
larger uncertainty of ±196(fcre/10keV) kms~ when 
the velocity is determined from combination of the 2.1 
and 1.4 mm bands. All three bands are used to deter- 
mine the peculiar velocity of A2163; we conservatively 
adopt the larger of the uncertainties determined from 
two bands. 

4. S-Z Data Analysis 

This section describes the reduction and analysis of 
the S-Z scans. The analysis of the A2163 2.1 mm data 
set is presented in detail in H97b. The initial reduc- 
tion of all the S-Z data is identical to that treatment. 
We outline the analysis method, with particular em- 
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702 


773 


151.5 
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420 


478 


-119.5 



average 261 550 627 



19.6 



Table 1: Residual thermal effect signal in the 1.4 mm band calibrated in terms of an "offset" in peculiar velocity. 
Results are shown for the three principal IC gas temperatures discussed in this paper. We also show the results 
using the non-relativistic S-Z spectra for 10 keV IC gas. 



phasis on the determination of accurate confidence 
intervals from multiparameter fits to the S-Z scans. 

4.1. Data Set 

Each of the two rows of the SuZIE array consists 
of three detectors; si, S2 and S3. Detector signals 
are differenced in pairs (in hardware) to form three 
difference signals; di2 — si — S2, ^23 = S2 — S3, and 
dsi — S3 — si. These differences correspond to angular 
chops of 2.3', 2.3', and 4.6' respectively. The two 2.3' 
differences are differenced to form a triple beam chop 
(TBC), for example, 



123 



di2 



-^23 



Si - 2S2 + S3 



(15) 



The data set consists of 4 difference signals: d^i, ti23, 
dQ4, and t456 corresponding to the 4.6' difference and 
TBC for each of the two rows of detectors. We refer 
to these four difference signals as dk, where k ranges 
from 1 to 4. We also ompute the average signal of the 
undifferenced detectors. 
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(16) 



The average single detector signal is used as a monitor 
of the absolute atmospheric emission and to remove 
any residual common mode response from the detec- 
tor differences. 

We clean the raw data of transients due to the 
interaction of cosmic rays with the detectors. Less 
than 5% of the data are identified as contaminated 
by cosmic rays. The raw data are binned into 3 s bins 



corresponding to 15 5 Hz samples or 0.75'cos((5) on 
the sky. Samples fiagged as bad due to cosmic rays 
arc left out when the bin averages are computed; bins 
with more than half of their samples fiagged are not 
used in the analysis. For each scan j, and bin i, the 
binned differential and average single detector signals 
are expressed as dkji and Sji. 

4.2. S-Z Surface Brightness Model 

The S-Z surface brightness profile of A2163 has 
a FWHM w 5', which is comparable to the largest 
beam-throw of SuZIE. In order to accurately deter- 
mine the peak surface brightness, we must simulate 
the observation of the extended source. Models for 
the surface brightness of the S-Z thermal and kine- 
matic components are constructed from the X-ray 
surface brightness determined density profile and the 
assumption of some thermal structure. We express 
the surface brightness morphology in terms of dimen- 
sionless form factors normalized to one at peak bright- 
ness: 
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(17) 
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To take differences in the array elements into account, 
SK,Ti are evaluated at the S-Z intensity weighted 
band center for each array element i. In the case of 
an isothermal IC gas, Snid, (j)) = SKi{d, <t>) for all i. 
The form factor for the desired S-Z component 
{SK,Ti) is convolved with beam-maps constructed 



from the voltage response of the detectors to scans 
over planets, Vpi. This forms a template for the re- 
sponse of each detector to the assumed surface bright- 
ness of the S-Z components, 



smK,Ti{0) = 

AIK,T^ 1 
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(19) 
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where flp is the solid angle subtended by the plan- 
etary calibrator. This expression assumes the plane- 
tary calibrator to be much smaller than the beams on 
the sky, true for Mars and Uranus. To calibrate the 
model, we determine the ratio of the relativistically 
correct S-Z brightness to the brightness of the planet. 



Ik.tj _ S Ii{v)lK,Tudv 
Ipt Jfi{u)Ip^du 



(20) 



where fi{v) is the measured spectral response and Ip^ 
is the intensity of the planetary calibrator. The S-Z 
brightness is calculated for y = 1 and Vr = lkms~^, 
so the results of the model fits are in convenient units. 
The models for the response of the single array 
elements are differenced to create models for the re- 
sponse of the detector differences. The ~ .05' resolu- 
tion differential source models are binned identically 
to the scan data to determine the model signal for 
each of the ~ .75' data bins. The binned differential 
models are designated as niki (RA) , where k is one of 
the four detector differences, i is the position (by bin 
number) in the scan, and RA is a offset of the model 
from the nominal X-ray determined position. 

4.3. Coadded Data 

We coadd the difference signals to create high sen- 
sitivity scans of the differential surface brightness 
as a function of RA. For each difference channel k, 
and scan j, we clean the data by removing the best 
fit linear baseline and residual common-mode signal 
[akjSji), 
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(21) 



The residual common-mode signal is the small signal 
in the detector differences that is correlated with the 
average single detector signal. Removing this signal 
improves the rejection of the differential signals to 
common mode atmospheric noise. In H97b, we show 
that removing it has no systematic effect on the fit 
results. 



The value of each coadded bin Xki, is given by the 
weighted sum of the value of this bin in each of Ng 
scans. 
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Each bin is weighted by the residual RMS of the scan, 
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where Nf, is the number of bins in a scan. The uncer- 
tainty in the value of each bin is determined from the 
weighted dispersion of the value of that bin about the 
mean determined from Ns scans, 
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The best fit peak Coniptonization (j/g) and isother- 
mal model position (RA) are found by minimizing the 
X^ of the fit to all four difference signals for each of 
the coadded 2.1 mm data sets. 
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(25) 
The linear baseline is unconstrained to make sure that 
our removal of a linear baseline in equation G^ did not 
effect the signal. For each value of RA, we re-bin the 
.05' resolution source model to determine the model 
signal for each of the ~ .75' data bins. 

If the source contributes significantly in a single 
scan, then a scan in which the noise and source con- 
spire to produce a low RMS and corresponding low 
signal will be weighted higher. Coadding the scans 
using equations |2l|-|2^ could introduce a bias in the 
amplitude of the coadded scans. In H97b, we elimi- 
nate this bias by recomputing the coadded bins with 
each bin weighted on the residual RMS with the av- 
erage source model removed. We do this for all the 
coadded data presented in this paper as well, although 
it is unlikely to be necessary for the extremely low sig- 
nal to noise 1.4 and 1.1 mm scans. 



4.4. Single scan fits 

The determination of accurate uncertainties for fits 
to the coadded scans is difficult due to the presence 
of atmospheric noise which is correlated between bins 
and detector differences. Ignoring this correlation 
leads to a significant overestimate of the significance 
of the fit results. In this section, we use the distri- 
bution of the single scan fit amplitudes to determine 
accurate uncertainties for model amplitudes. These 
results are then used to correct the error in the coad- 
ded fit uncertainties due to correlated sky noise. 

The noise in scans taken at different times is uncor- 
related. Each scan can be treated as an independent 
measurement of the source amplitude. We determine 
the mean single-scan peak Comptonization for a given 
observation by averaging the values yoj for each of the 
Ns scans, 
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where ayj is the change in j/oj corresponding to Ax^ = 
1 in the model fit. The weighted dispersion of the 
scan amplitudes about the mean is used to estimate 
the uncertainty in the determination of the mean, 



\ 



E 




(Ns- 


Ns 



(27) 



In H97b, we plot the distribution of single scan fit 
amplitudes for all of the 2.1 mm scans across A2163. 
As in that example, the distribution of scan fit am- 
plitudes for all the data sets considered here are well 
approximated by a Gaussian. The width of the dis- 
tribution is a function of the amount of sky noise and 
is correlated with the RMS of the scans. Using the 
weighting scheme described in this section, we narrow 
the distribution slightly by decreasing the weight of 
the scans with large sky noise. 

4.5. Confidence intervals 

From the single-scan fits we have an accurate mea- 
surement of the uncertainty in the source amplitude 
for each data set, <Ty. We compare this to the uncer- 
tainty (cTyca) determined from the fits to the coadded 
data (equation p3; Ax^ = I). We compute a scaling 



factor to correct for the underestimation of the bin 
uncertainties for each data set, m. 
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Depending on the details of the atmosphere at the 
time of the measurement, 7„i = 1.4 — 2.0. Each coad- 
ded scan is individually scaled so that the uncertain- 
ties in the coadded scan fit amplitudes are equal to the 
single-scan fit uncertainties. Fits to the coadded data 
(for similar model shapes and positions) will then re- 
sult in accurate uncertainties for the fit parameters. 
Essentially, we have reduced the number of degrees 
of freedom to account for the correlated noise. The 
resulting confidence intervals are a factor of 1.4 to 2.0 
larger than those found from the uncorrected data. 
All confidence intervals in this paper are computed 
using coadded data that has been corrected for cor- 
related sky noise in this way. 

We determine confidence intervals for the source 
position, peculiar velocity, peak Comptonization and 
isothermal P model parameters using a maximum 
likelihood estimator. 
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where Cfc = Xki - yo mTkt{.RA, 0^ (3) 
— Vr mKki{RA, 6, (3) — Qk — i bk is the coadded scan 
with the kinematic and thermal component source 
models removed. The coadded scans corresponding 
to different RAOs and frequency bands are treated as 
independent; the likelihood for simultaneous fits of k 
scans for each of m uncorrelated data sets is given by. 



^(o=nn^'"fc(c) 



(30) 



To make the problem computationally tractable we 
determine the likelihood for only two interesting pa- 
rameters at a time. The likelihoods are determined 
for a necessarily large grid in the parameter space 
with resolution: ARA — .05', Avr — lOkms^ , and 
Aj/o = 2 X 10^^. Assuming no peculiar velocity, fits 
to the X-ray derived isothermal thermal component 
model determine the source position and amplitude. 
The resulting position is insensitive to the values of 
the peculiar velocity and /3 model parameters. For 
the rest of the analysis, we fix the source position to 
its best fit value. The isothermal /3 model parameters 



are determined with the source amplitude free to vary. 
We do this to check that the assumed model is con- 
sistent with the S-Z data. Finally, we fix the position 
and (3 parameters for the the thermal and kinematic 
component models and determine the likelihood for 
Vr and j/o- 

The likelihood grids are converted to confidence re- 
gions and intervals in the interesting parameters. In- 
voking Bayes' theorem and assuming a uniform prior, 
the probability that the two fit parameters (y, z) fall 
within a region R is given by 



P{R)^ 



I^oo I^oo ^(y^ ^) '^y^^ 



(31) 



The confidence interval corresponding to a probability 
Pq is given by the region Rq such that P{Ro) = Pq 
and L[{y, z) e Rq] > L[{y, z) ^ Rq]. We determine 
the likelihood in a single parameter by marginalizing 
the likelihood over the uninteresting parameters, 
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(32) 



The probability that the parameter y falls within the 
interval I — [zi, Z2] is given by, 



P{I) = 



Ii L{z) dz 

r^L{z)dz 



(33) 



The confidence interval corresponding to a probability 
Pq is given by Iq — [zi, Z2] such that P{Iq) — Pq and 
i(-z G 4) > i(2 ^ -^- AH the 68% confidence levels 
quoted in Sections ^ and |6j are calculated this way 
unless otherwise stated. 

4.6. Baseline 

We have carefully designed our instrument and 
scan strategy to eliminate systematic errors. How- 
ever, it is possible that the instrument could introduce 
some systematic signal similar to the expected cluster 
signal. To test for the presence of this type of instru- 
mental baseline, we gather "blank sky" scans across 
regions of sky free of known sources. These data are 
analyzed exactly as the source data; the presence of 
an instrumental baseline would result in the determi- 
nation of a non-zero signal for these regions. 

In April 1994, we accumulated ^ 15 hours of 
2.1mm data on two patches of sky free of known 
sources. The analysis of these data are described 
in Holzapfel et al. (1997b). Fitting these data to 



the isothermal thermal effect model for A2163 with 
the position in the scan fixed at the best fit po- 
sition from the source observations, we determine 
2/0 = 1-55 ± 2.13 X 10"^ In May 1993, we accumu- 
lated ^ 7.2 hours of 1.1 mm data on three patches of 
sky free of known sources. Fitting these data to the 
isothermal thermal effect model for A2163, we deter- 
mine 2/0 = -.21 ±2.17 X 10"**. In May 1994, we accu- 
mulated ~ 16.3 hours of 1.4 mm data on two patches 
of sky free of known sources. Fitting these data to the 
kinematic component source model for A2163 and as- 
suming T — .015, we find Vr — -1-760 ± 1050 kms~ . 
Therefore, we detect no significant baseline in any of 
the three spectral bands. 

We could subtract these scans from the source 
scans to remove a possible instrumental baseline. Un- 
fortunately, due to shorter integrations and poorer 
weather, the "blank sky" data in the 1.4 and 1.1mm 
bands is less sensitive than the source data; the sub- 
traction of these scans would significantly increase the 
uncertainty in the fit results. Another problem with 
this procedure is that the baseline and source data, al- 
though gathered at similar azimuth and zenith angle, 
are gathered at different times. Any baseline signal, 
because it is not significant in a single scan, must be 
correlated in time between several scans. However, 
this does not guarantee an instrumental baseline will 
be constant over the course of the night. 

We have devised an observation strategy which al- 
lows us to check for, and remove, any time correlated 
baseline. Differencing scans taken adjacent in time 
at two different RAOs allows the subtraction of any 
common baseline while retaining most of the expected 
source signal. Each pair of scans are differenced by 
subtracting their raw 5 Hz sampled time streams. The 
data are then analyzed exactly as described in Sec- 
tions |4]-|4.5. 



We construct a model for the RAO differenced sur- 
face brightness by subtracting two models with the 
appropriate RAO difference between them. We chose 
the RAO difference (6 ) to be large enough so that 
(for most clusters) the two models have little over- 
lap in time and therefore will not subtract much sig- 
nal when differenced. The RAO differenced model 
is fit to the RAO differenced data to determine the 
source amplitude and position. In all cases the results 
agree, within the quoted statistical error, with those 
found from fits to both RAOs. The lack of any signif- 
icant change when the scans are differenced indicates 
that there is no significant instrumental baseline in 



the data. We use the difference between the results 
of the RA differenced and standard analysis as an es- 
timate of the uncertainty due to the possibility of an 
instrumental baseline. 

5. A2163 

5.1. X-ray Surface Brightness 

Abell 2163 is a distant (z = .201) X-ray lumi- 
nous cluster of galaxies with the distinction of having 
the highest X-ray temperature of any cluster yet ob- 
served. A2163 was the target of pointed observations 
by the GINGA, ROSAT and ASCA satellites. A de- 
tailed analysis of each of these observations has been 



^ 



ublished (|Arnaud et al 1992| ; |Elbaz et al. 1994 



arkcvitch et al. 1996|) . The ROSAT/PSPC was 



used to determine the spatial dependence of the X-ray 
surface brightness. The peak surface brightness was 
found at 16^15'"46^ -06° 09' 02" (J2000). The ra- 
dial profile was determined by summing annuli about 
this position, with significant emission detected up 
to 18' from the position of the peak. Assuming the 
IC gas to be isothermal, the surface brightness profile 
was fit with the combination of an isothermal (3 model 
and a constant background. The density profiles are 
assumed to have the form of a spherically symmetric 
/? model. 
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The best fit model parameters were found to be /3 = 
.61 6tmq and Oc = 1.20 ± 0.05' at 68.3% confidence 



( lElbaz et al. 1994]) 



5.2. Thermal Structure 

In H97b, the thermal structure of the IC gas in 
A2163 is discussed in detail, here we summarize only 
the essential points. A combined analysis of the 
ASCA/GIS+SIS and GINGA spectral data gives tem- 
peratures: 13.3, 13.3, and 3.8 keV in three rad ial bins 
of — 3, 3 — 6, an d 6 — 13 X-ray core radii ( Marke- 
dtch et al. 1996). Wc use a simple two parameter 



model to describe the measured temperature profile 
of A2163. The model consists of an isothermal central 
region extending to 9iso beyond which the tempera- 
ture decreases according to a polytropic model with 



index 7, 
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The central temperature is Teo ~ 13.3^^gkeV, where 
the error bars enclose the 68% confidence intervals for 
both the GINGA and ASCA results, each analyzed in- 
dependently under the assumption of the ASCA ther- 
mal structure (II97b). We adopt thermal model pa- 
rameters Oiso — 4.0 9c and 7 = 2.0 and reanalyze the 
X-ray data to determine consistent density model pa- 
rameters /3 = .64 and Oiso = 1-26. These results are 
used to create a model for the thermal component sur- 
face brightness that includes the thermal structure. 

In the case when the IC gas is assumed to be 
isothermal, we adopt a temperature which is the 
emission-weighted average temperature of the above 
results. For an isothermal IC gas, Tg « 12.4^^7keV, 
where the error bars enclose the 68% confidence in- 
tervals for both the GINGA and ASCA results, each 
independently analyzed assuming an isothermal IC 
gas. The results in this section are for fits to the 
isothermal source model unless otherwise stated. 

5.3. 2.1 mm Observations 

We have observed A2163 at 2.1 mm in two periods; 
April 23-26 1993 for a total of 16 hours and April 
4, 9, 10, and 11 for 8 hours. The analysis of these 
data to determine the Hubble constant is presented 
in Holzapfel et al. (1997b). 

The data corresponding to each observation run 
and RAO are combined to make coadded scans. In 
Figure g, we plot the April 1994 coadded 4.6' and 
TBC scans over A2163 with the best fit isothermal 
model overlaid. Assuming the peculiar velocity to be 
zero, the coadded data scans are simultaneously fit to 
the corresponding isothermal source models with the 
peak S-Z surface brightness position and peak Comp- 
tonization free to vary. In Figure ^, we show the 68.3 
and 95.4% confidence regions for the source position 
and amplitude [RA and r/o)- The 68.3% confidence 
limits for the fit parameters considered individually 
are /\RA = .35 ± .14' and yo = 3.73 ± 0.35 x lO"'*. 
Including the uncertainties in the ROSAT/PSPC and 
SuZIE pointing, the positions of the peak X-ray and 



S-Z surface brightnesses are consistent (Holzapfel ex 
al . 1997b| ). 

Fixing the position to the best fit value, we fit 
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Fig. 2. — Coadded data for all April 1994 scans across 
A2163 at 2.1 ,mm. The coadded scans from the two 
RAOs have been offset and added to create a single 
coadded scan for each difference. Both the 4.6' dif- 
ference (upper panel) and TBC (lower panel) data 
are shown for the row passing directly over the peak 
X-ray surface brightness. The heavy line is the best 
fit isothermal model, and the lighter lines represent 
the 1(7 errors. The scans are calibrated in terms of 
a Rayleigh- Jeans temperature difference that com- 
pletely fills one of the beams. 

isothermal /? models for the surface brightness dis- 
tribution to the coadded scans. In Figure |j, we plot 
approximate confidence regions in 9c and j3. Techni- 
cally, the data lack sufficient sensitivity to determine 
proper confidence regions in position and amplitude; 
we approximate the 68.3 and 95.4% confidence regions 
as corresponding to Ax^ — 2.3 and 6.17 respectively. 
The results of the S-Z fits are consistent with those of 
the X-ray derived density parameters and an isother- 
mal IC gas. In H97b, we found that the S-Z data 
did not have sufficient sensitivity to discriminate be- 
tween an isothermal IC gas and the thermal structure 
measured with ASCA. 

5.4. 1.1 mm Observations 

We observed A2163 on May 1-4 1993 for a total of 
19.7 hours. Scans were 30' long with one row pass- 
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Fig. 3. — Isothermal thermal component model model 
fit to the coadded 2.1 mm data for A2163. Assuming 
no peculiar velocity, the contours enclose the approx- 
imate 68.3% and 95.4% confidence regions in peak 
Comptonization and RA offset from the X-ray cen- 
ter. 

ing over the X-ray center and the other 2.2' to the 
south. The RAO of the X-ray center from the start 
of the scan was alternated between 12.5' and 18.5' in 
sequential scans. The 1.1mm observations are domi- 
nated by sky noise and the quality of the data varies 
greatly with atmospheric conditions. 

In Figure S, we show the coadded 4.6' and TBC 
scans over A2163 with the best fit isothermal model 
overlaid. Assuming the peculiar velocity to be zero, 
the coadded data scans are simultaneously fit to the 
corresponding isothermal source model with the po- 
sition and peak Comptonization free to vary. In Fig- 
ure p, we show the approximate confidence regions in 
RA and j/o computed for the fits of the isothermal 
source model to the 1.1mm coadded data. The ap- 
proximate single-parameter 68.3% (Ax^ = 1.0) con- 
fidence intervals are /S.RA = .17 ± .59' and j/o = 
4.05 ± 1.47 X 10""'. Both the amphtude and posi- 
tion of the peak S-Z surface brightness are consistent 
with the results of the 2.1 mm analysis. 

5.5. 1.4 mm Observations 

We observed A2163 in the 1.4 mm band on the 
nights of May 5, 6, 7, and 8 1994 for a total of 
~ 12.3 hours. Scans were 30' long with one row pass- 
ing over the X-ray center and the other 2.2' to the 
south. The RAO of the X-ray center from the start 
of the scan was alternated between 12' and 18' in se- 
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Fig. 4. — Isothermal density model parameter fits to 
the coadded data scans across A2163. Contours for 
6c and (3 are shown at approximate ly 68% and 95% 
confidence. The "X" marks the best fit isothermal 
model parameter valu es from the X-ray analysis. 

quential scans. In Figured, we show the coadded 4.6' 
and TBC scans. Fitting the data with a model for the 
kinematic component, we determine a marginally sig- 
nificant signal, the bulk of which is due to the residual 
thermal component signal in the 1.4 mm band. 

5.6. Simultaneous Fits 

5.6.1. 2 Band Fits 

We simultaneously fit the 2.1 and 1.1mm coad- 
ded scans with the S-Z thermal and kinematic com- 
ponent isothermal models generated for each of the 
two bands. The models are fixed to the best fit po- 
sition from the 2.1mm scans. We determine Vr = 
-l'20tlllo kms"^ and yo = 3.92 ±.59 x IQ-^ at 68.3% 
confidence. 

We also simultaneously fit the 2.1 and 1.4 mm 
coadded scans with the S-Z thermal and kinematic 
component isothermal models generated for each of 
the two bands. The models are fixed to the best 
fit position from the 2.1mm scans. We determine 



Vr = +990+1190 kms" 
68.3% confidence. 



and yo = 3.45 ± .63 x lO"* at 



5.6.2. 2.1 mm, l.Am,m, and 1.1 mm Bands 

We simultaneously fit isothermal (Tg = 12.4 keV) 
models for the kinematic and thermal components of 
the S-Z effect to the 2.1, 1.4 and 1.1mm coadded 
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Fig. 5. — Coadded data for all May 1993 scans across 
A2163 at 1.1mm. The scans for the two RAOs have 
been offset and added to create a single coadded scan 
for each difference. Both the 4.6' difference (upper 
panel) and TBC (lower panel) shown for the row pass- 
ing directly over the peak X-ray surface brightness. 
The heavy line is the best fit isothermal model, and 
the lighter lines represent the la errors. The scans are 
calibrated in terms of a Rayleigh-Jeans temperature 
difference that completely fills one of the beams. 

scans. The model positions are fixed to the best fit po- 
sitions from the 2.1 mm fits; the model amplitudes, j/o 
and Vr^ are left as free parameters. We plot the 68.3% 
and 95.4% confidence regions for j/q and v,. in Figure pi 
The 68.3% confidence intervals for the fit parameters 
considered individually are Vr = -I-490I730 kms~^ and 
2/0 = 3.62 ± .48 X 10"''. 

We have repeated this analysis for the range of 
allowed density models; the results are unchanged 
within the precision of the grid in Vr (10 kms~ ) . This 
analysis is repeated using the RAO differenced scans; 
we compare the results of the two analysis methods 
in Table g. The sensitivity in v,. is much worse for the 
RAO differenced data. This is due to the large extent 
of the A2163 surface brightness. A larger RAO dif- 
ference would reduce the amount of subtracted signal 
and increase the sensitivity of the RAO differenced 
data. 
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Peculiar Velocities: Standard and RAO Differenced 



Cluster Thermal Structure fcTgo [keV] analysis yo x 10"* Vr [kms ^] 



A2163 


Isothermal 


A2163 


Isothermal 


A2163 


ASCA 


A2163 


ASCA 


A1689 


Isothermal 


A1689 


Isothermal 



12.4 


Both RAOs 


3.62 ± 


.49 


+49Gt?J° 


12.4 


Differenced 


3.44 ± 


.72 


+ 1390i2350 


13.3 


Both RAOs 


2.96 ± 


.39 


+560t^^g 


13.3 


Differenced 


2.34 ± 


.57 


+1460l?g° 


8.2 


Both RAOs 


3.43 ± 


.59 


+17011?;; 


8.2 


Differenced 


3.00 ± 


.73 


+310l^r 



Table 2: Comparison of peculiar velocities determined with the standard analysis and fits to the RAO differenced 
scans. All results are at 68% confidence. The differences between the results represent our best limits on possible 
contribution due to an instrumental baseline. The sensitivity is considerably worse for the RAO differenced scans, 
especially for A2163, which is more extended than A1689 and has more signal subtracted in the differencing. 
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Fig. 6. — Isothermal thermal component model model 
fit to the coadded 1.1 mm data for A2163. Assuming 
no peculiar velocity, the contours enclose the approx- 
imate 68.3% and 95.4% confidence regions in peak 
Comptonization and RA offset from the X-ray center 
(Note: scale is different than Figure g). 



We have repeated the determination of v.^ assuming 
two other values for the isothermal IC gas tempera- 
ture which roughly span the 68.3% confidence interval 
for Tg. The peculiar velocities determined assuming 
Te = 10 and 15 keV are listed in Table | The effects 
of the implied change in optical depth and realtivis- 



tic corrects cancel. The contribution of the uncer- 
tainty in the IC gas temperature to the uncertainty 
in peculiar velocity {Avr < ±20kms~^) is negligible 
compared to the statistical uncertainty in the mea- 
surement. 

5.7. 2.1, 1.4, and 1.1mm Fits With Thermal 
Structure 

ASCA observations of A2163 suggest that the IC 
gas is not isothermal. We have repeated the simul- 
taneous fits to the A2163 S-Z data for the case in 
which the cluster is assumed to have the ASCA de- 
termined thermal structure (Section 5.2.). Using the 



hybrid IC gas temperature model and the consis- 
tent density model, we create models for the surface 
brightness of the thermal and kinematic components 
of the S-Z effect. For a single free parameter, we find 
Vr = -t-560l^[;^kms"^ and yo = 3.00 ± .40 x 10"^ at 
68.3% confidence. These values are compared with 
the isothermal results in Table 0. We have repeated 
this analysis with the RAO differenced scans; We 
compare the results of the two analysis methods in 
Table Although there is considerable uncertainty in 
the thermal structure of the cluster, the uncertainty 
it contributes to the determination of the peculiar ve- 
locity is negligible in comparison to the statistical un- 
certainty. 
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Peculiar Velocities: IC gas Temperature 



Cluster 


Thermal Structure 


kTeo [kcV] 


yo X 10^ 


Vr [kms ^] 


A2163 


Isothermal 


10.0 


3. 50 ±.47 


+510t^« 


A2163 


Isothermal 


12.4 


3.62 ±.49 


+490j:?j° 


A2163 


Isothermal 


15.0 


3. 70 ±.50 


+mtl°l° 


A2163 


Hybrid(ASCA) 


13.3 


2.96 ±.39 


+560t^g° 


A1689 


Isothermal 


5.0 


3.31 ± .52 


+190i^™ 


A1689 


Isothermal 


8.2 


3.43 ±.59 


+170^^7^ 


A1689 


Isothermal 


10.0 


3. 51 ±.55 


+130t^?° 



Table 3: Peculiar velocities assuming several values of the IC gas temperature for each cluster. For A2163, we 
also show the results for the Hybrid(ASCA) thermal structure. Results are for fits to the coadded data at 68% 
confidence. 



6. A1689 

6.1. X-ray Results 

AbeU 1689 is a distant {z = .181) and X- 
ray luminous cluster of galaxies. It has been ob- 
served with the ROSAT, GINGA, and ASCA satel- 
lites. ROSAT/PSPC observations determined the 
peak of the X-ray surface brightness to be centered 
at 13'^11'"29.1^ -01° 20' 41'^ (J2000) with sign ificant 
emission detected to 8.5' ( Daines et al. 1995| ). The 
X-ray surface brightness contours are nearly circu- 
lar with ellipticity < 0.1. Assuming the density of 
the IC gas to be described by an isothermal beta 
model (equation M), Daines et al. (1995) determine 
0c = 1.13 ± 0.12' and f3 = 0.78 ± 0.03 at 90% confi- 
dence. From the observed radial velocities of the clus- 
ter member galaxies, it appears there are two to three 
substructures superimposed along the line of sight. 
While this may complicate the use of this cluster to 
determine the Hubble constant, it has no effect on the 
determination of the peculiar velocity. 

The presence of a cooling flow (M « 500MQyr^^) 
within the central core of the cluster complicates the 
X-ray spectral analysis. Daines et al. (1995) use the 
combined analysis of GINGA and ROSAT spectra to 
determine T^ ^ 9 keV outside the cooling flow. Anal- 
ysis of ASCA spectra has also been used to determine 
the temperature of the IC gas. Assuming the gas to 
be isothermal outside of the central core region the 



t emperature is det ermined to be Tg 
( iBautz et al. 1994| ). 



8.2 ± 1.0 keV 



6.2. 2.1 mm Observations 

We observed A1689 in the 2.1mm band on the 
nights of April 6, 7, and 9 1994 for a total of ~ 
8.5 hours. Scans were 30' long with one row passing 
over the X-ray center and the other 2.2' to the south. 
The RAO of the X-ray center from the start of the 
scan was alternated between 12' and 18' in sequential 
scans. 

In Figure ^ we show the coadded 4.6' and TBC 
scans over A1689 in Figure 0. Assuming no peculiar 
velocity, we fit the isothermal model for the S-Z ther- 
mal component to the coadded scans. In Figure |lO|, 
we plot the 68.3 and 95.4% confidence regions for yo 
and RA. The 68.3% confidence intervals for the posi- 
tion and peak Comptonization are ARA = —.01 ±.11' 
and yo = 3. 55 ±.29 x 10^^. The positions of the X-ray 
and S-Z surface brightnesses are coincident. 

Combining the peak Comptonization with the X- 
ray determined IC temperature, we determine the op- 
tical depth, r ~ 0.022. This extremely high value of 
optical depth is exceptional and suggests the exten- 
sion of this cluster along the line of sight. The sensi- 
tivity of our measurements to cluster peculiar velocity 
scale linearly with r making this an ideal cluster for 
the determination of the peculiar velocity. 

Fixing the position to the best fit value, we fit 
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Fig. 7. — Coadded data for all May 1994 scans across 
A2163 at 1.4 mm. The coaded scans from the two 
RAOs have been offset and added to create a single 
coadded scan for each difference. Both the 4.6' differ- 
ence (upper panel) and TBC (lower panel) are shown 
for the row passing directly over the peak X-ray sur- 
face brightness. The heavy line is the best fit isother- 
mal model, and the lighter lines represent the la er- 
rors. The scans are calibrated in terms of a Rayleigh- 
Jeans temperature difference that completely fills one 
of the beams. 

isothermal f3 models for the surface brightness dis- 
tribution to the coadded scans. In Figure ^ We plot 
approximate confidence regions for 6c and (3. The ob- 
served S-Z surface brightness is consistent with the 
best fit isothermal model parameters from the X-ray 
analysis and the assumption of an isothermal IC gas. 

6.3. 1.4 mm Observations 

We observed A1689 in the 1.4 mm band on the 
nights of May 6-9 1994 for a total of ~ 14 hours. The 
scans were 30' long with one row passing over the X- 
ray center and the other 2.2' to the south. The RAO 
of the X-ray center from the start of the scan was al- 
ternated between 12' and 18' in sequential scans. In 
Figure |l^, we show the coadded 4.6' and TBC scans 
over A1689. Fitting the data to a model for the kine- 
matic component of the S-Z effect, we determine a 
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Fig. 8. — Results of fitting the 2.1, 1.4, and 1.1mm 
coadded data for A2163 to isothermal (Tg = 12.4 keV) 
source models for the two components of the S-Z ef- 
fect. The 68.3% and 95.4% confidence regions in peak 
Comptonization and radial component of the peculiar 
velocity are shown. 



marginally significant signal, the bulk of which is due 
to the residual S-Z thermal component. 

6.4. Simultaneous Fits 

We simultaneously fit isothermal (Te = 8.2 keV) 
models for the kinematic and thermal components of 
the S-Z effect to the 2.1 and 1.4 mm coadded scans. 
The model positions are fixed to the best fit positions 
form the 2.1mm fits; the model amplitudes, j/o and 
Vr, are left as free parameters. In Figure O, we plot 
the 68.3% and 95.4% confidence regions for yo and Vr- 
The 68.3% confidence intervals for the fit parameters 
considered individually are yo = 3. 43 ±.59 x 10^* and 
Vr = -l-170l57Qkms~ . The results are unchanged 
for the allowed range of density model parameters. 
We have repeated this analysis with the RAO differ- 
enced scans; the results of the two analysis methods 
are compared in Table g. 

We have determined yo and Vr assuming two other 
values for the isothermal IC gas temperature which 
exceed the range of temperatures allowed at 68.3% 
confidence. The peculiar velocities determined as- 
suming Te — 5 and lOkeV are listed in Table |[ 
The contribution of the uncertainty in the temper- 
ature to the uncertainty in peculiar velocity (Aw^ < 
±50 kms^ ) is much smaller than the statistical un- 
certainty in the measurement. 
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Fig. 9. — Coadded data for all April 1994 scans across 
A1689 at 2.1 mm. The coaded scans from the two 
RAOs have been offset and added to create a sin- 
gle coadded scan for each difference. Both the 4.6' 
difference (upper panel) and TBC (lower panel) data 



arel shown tor the row passing directly over the peak 



X-ray surtace brightness. I'he heavy line is the best 
fit isothermal model, and the lighter lines represent 
the la errors. The scans are calibrated in terms of 
a Rayleigh- Jeans temperature difference that com- 
pletely fills one of the beams. 

7. Astrophysical Confusion 

Although astrophysical confusion from randomly 
distributed sources is expected to be small at mm 



wavelengths (Fischer and Lange 1993), the possibil- 
ity of confusion contributing a considerable system- 
atic error to the determination of Vr in an individual 
cluster must be considered. In the section, we exam- 
ine the contribution of confusion, including primary 
anisotropics of the CMB, to the uncertainty in pecu- 
liar velocity for the two clusters considered. 

7.1. Radio Source Confusion 

7.1.1. A2163 

A VLA search towards A2163 shows evidence of a 
radio source 0.8' west of the cluster center with an 
inverted spectrum (Herbig fc Birkinshaw 199^). For 
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Fig. 10. — Isothermal thermal component model fit 
to the coadded 2.1 mm data for A1689. Assuming no 
peculiar velocity, the contours enclose the 68.3% and 
95.4% confidence regions in peak Comptonization and 
RA offset from the X-ray center. 

this source, the flux rises from 1 to 3 mJy between 
6 and 2 cm suggesting a flux as large as 30 mJy at 
2.1mm. Recently, this cluster has been imaged with 



the OVRO interferometer operating at 30 GHz (Carl- 



Strom et al. 1996). There is only one significant un- 
resolved source found in the field. It is centered at 
!&", IS"", 43.7^ -06°, 08', 45" (J2000), close to the 
VLA source position, and has a flux of 1.4 mJy. Due 
to the steeply falling spectra from 2 cm to 1 cm, the 
source flux is likely to be less than 1 mJy in any of 
the SuZIE bands. However, this sou rce is reported to 
be time variable ( Birkinshaw 1996 ). Because of the 
elapsed time between measurements, the spectral in- 
dex and the inferred fluxes in the SuZIE bands are 
questionable. The central region of the cluster has 
been imaged at 3.3 mm with the IRAM interferome- 
ter. Fischer & Radford report an upper limit of 5 mJy 
(2tT in a 20" (5, 10" RA beam) on point source emis- 
sion within 1' of the X-ray center. More sensitive 
observations with mm- wavelength interferometers are 
needed to eliminate the possibility of confusion due to 
this particular radio source. At present, we use the 
measured spectral index to infer that the fiux in the 
SuZIE bands is likely to be less than 1 mJy. 

A2163 hosts the brightest radio halo yet discovered 
(Herbig & Birkinshaw 1995). From measurements at 
1.5 and 4.9 GHz, the integrated flux from the radio 
halo is estimated to be less than 1 mJy in the SuZIE 
2.1 mm band (Herbig 1995|). Taking the measured op- 
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Fig. 11. — Isothermal density model parameter fits to 
the coadded data scans across A1689. Contours for 
6c and (3 are shown at approximately 68% and 95% 
confidence. The "X" marks the best fit isothermal 
model parameter values from the X-ray analysis. 

tical depth of A2163 and beam dilution into account, 
a flux of 1 mjy in any of the three observing bands 
corresponds to an error of ~ 120 kms~ in the derived 
peculiar velocity for A2163. The total uncertainty in 
the peculiar velocity of A2163 due to radio confusion 
is estimated to be < ±170kms^^. 

7.1.2. A1689 

Radio emission from A1689 has been studied in 
detail by Slee et al. (1994). They find two bright 
sources near the cluster center which they label 
#13 and #15a/b. Source #13, at central position 
13^^, 11'", 

30.2''; -01°, 20', 29" (J2000), has a flux of 10.4mJy 
and 2.75 mJy at frequencies of 1.5 GHz and 4.9 GHz. 
Source #15a/b, at central position 13'', 11'°, 31.4"; 
-01°, 19', 33" (J2000), has a flux of 41.1 mJy and 
10.7 mJy at frequencies of 1.5 GHz and 4.9 GHz. Us- 
ing these results, they determine the spectral indices 
to be a = —1.12 and a = —1.14 for sources #13 and 
#15a/b respectively. 

Recently, sensitive maps of A1689 have been made 
with the OVRO interferometer operating at 30 GHz 
(Carlstrom et al. 1996). They find two signifi- 




cant point sources in the field with approximate po- 
sitions 13'', 11"', 30. 3^ -01°, 20', 26" (J2000) and 
13'', ll"", 31.0^ -01°, 19', 36" (J2000) and fluxes 1.3 
and 0.4 mJy. Within the limits of the accuracy of 
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Fig. 12. — Coadded data for all April 1994 scans 
across A1689 at 1.4 mm. The coaded scans from the 
two RAOs have been offset and added to create a sin- 
gle coadded scan for each difference. Both the 4.6' 
difference (upper panel) and TBC (lower panel) data 
are shown for the row passing directly over the peak 
X-ray surface brightness. The heavy line is the best 
fit isothermal model, and the lighter lines represent 
the la errors. The scans are calibrated in terms of 
a Rayleigh- Jeans temperature difference that com- 
pletely fills one of the beams. 

the measurement (~ ±5"), these positions are coinci- 
dent with the previously identified radio sources. The 
fluxes are used to determine the spectral indices from 
4.9 to 30.0 GHz. We find a = -1.06 and a = -1.14 
for sources #13 and #15a/b, demonstrating that they 
continue their steep decrease in flux to 30 GHz. 

Using the measured spectral indices, we estimate 
the flux in all three SuZIE bands to be less than 
0.3 mJy from the brighter of the two sources. This is 
especially conservative in that the brighter of the two 
sources is « 1' north of the scan path, and therefore 
is never inside the FWHM of any of the array beams. 
The smaller beam dilution and higher optical depth of 
A1689 means that a 1 mJy point source in each of the 
three SuZIE bands produces only « 63kms~ error. 
We conclude that radio source confusion contributes 
< ±19 kms~ to the determination of the peculiar ve- 
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Fig. 13.— Results of isothermal (Te = 8.2 keV) model 
fit to correlation corrected corrected 2.1 and 1.4 mm 
coadded scans over A1689. The 68.3% and 95.4% 
confidence regions in the peak Comptonization and 
radial component of the peculiar velocity are shown. 

locity in A1689. 

7.2. IR Cirrus 

We have used IRAS 100 /im sky maps to place 
limits on confusion due to IR cirrus emission associ- 
ated with our galaxy. The images are obtained from 
the Infrared Sky Survey Atlas (ISSA) made available 
over the internet through the Infrared Processing and 
Analysis Center (IPAC). The images have had zodi- 
acal light subtracted and been "de-striped". Each 
image is 2° x 2° centered on the peak X-ray Surface 
Brightness for each cluster. They are composed of 
1.5' square pixels, although the actual resolution of 
IRAS at 100 /im is approximately 4' x 5'. 

The line of sight to A2163 passes near the edge of 
the galactic plane. In our scan region, the average 
100 Atm flux is 17.2MJysr-i. The RMS variation of 
the flux in the 1.5' square pixels across the scan is 
0.78MJysr~^. When the map is convolved with the 
SuZIE differential beam patterns, we find that the 
100 /im differential flux at any point in the scan is 
< ±126mJy/beam. After the removal of a gradient 
in brightness, this result holds for fits to both the 4.6' 
difference and the TBC. In the analysis, we have not 
fixed the model position to a particular point in the 
scan, doing this will result in a lower limit. 

The line of sight to A1689 passes though a region of 
low IR cirrus emission. In our scan region, the average 



100 fj,m flux is 2.9 MJy sr^^ . The RMS variation of the 
flux in the image 1.5' square pixels is 0.16MJysr~-'^. 
When the map is convolved with the SuZIE beam 
patterns, we find that the 100 /zm differential flux at 
any point in the scan is < ±76 mJy/beam. This result 
appears to be dominated by residual striping in the 
map. A more sophisticated analysis may result in 
significantly lower limits. 

We scale the maps to our observing frequencies 
by assuming the emission in the 100 /im map is due 
entirely to 23 K dust with an emissivity, e ex v~^'^^ 
( [Wright et al. 19911 ). In Table |, we list the upper 
limits on the differential fiux from IR cirrus scaled in 
the SuZIE bands. The fiux limits are used to calcu- 
late the maximum confusion in peculiar velocity due 
to anisotropic dust. 



Upper Limits for IR Cirrus Confusion 



Cluster 



Wavelength 
2.1 1.4 1.1 



0.17 


0.85 


1.50 


20 


102 


180 


0.10 


0.51 


0.93 


6 


32 


59 



A2163 Flux [mJy/beam] 
Avpec [kms^ ] 

A1689 Flux [mJy/beam] 
Avpec [kms""^] 



Table 4: Upper limits on peculiar velocity confusion 
due to IR cirrus in the SuZIE bands. A1689 was ob- 
served only in the 2.1 and 1.4 mm bands. 

Cluster member galaxies with slowly rising spectra 
could make a significant contribution to the uncer- 
tainty in peculiar velocity. Observations with exist- 
ing mm- wavelength interferometers could reduce this 
uncertainty for a r = .01 cluster to < lOOkms^ in a 
single 1/2 day observation. 

7.3. Primary Anisotropies 

It is also possible for the measurement of the S- 
Z effect to be confused by the presence of primary 
CMB anisotropies. The spectrum of these distortions 
is identical to that introduced by the S-Z kinematic 
effect and, therefore, especially serious for the de- 
termination of cluster peculiar velocities. The con- 
fusion in Vr, due to primary anisotropies of am- 
plitude AT/T, is Ai;^ = 100 (.01/r) [(AT/r)/3 x 
10~^] kms^^. Haehnelt and Tegmark (1996) have esti- 
mated the confusion limits from primary anisotropies 
to the determination of cluster peculiar velocities. For 
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Q = 1 (CDM) models with iltaryon = -01 — .1, cluster 
optical depth r = 0.015, and the SuZIE beam size, 
they find |Aupec| < 300 kms"^. This contribution to 
uncertainty as well as those due to foreground sources 
are listed for both clusters in Table 0. 

Peculiar Velocities and Uncertainties 



Uncertainty 



A2163 A1689 



Statistical 

Calibration (Flux) 

Calibration (Spectral) 

Baseline 

IC Gas Temperature 

Thermal Structure 

Radio Confusion 

IR Cirrus 

Primary Anisotropics 



400+910 

+ 174 

-174 

+240 

-240 

+900 

-0 

+20 

-10 

+70 

-0 

+ 170 

-170 

+ 180 

-180 

+300 

-300 



170 



-760 

-570 

+5 

-5 

+ 160 

-160 

+ 140 

-0 

+20 

-40 



+ 19 
-19 
+32 
-32 
+200 
-200 



Total 



490: 



-1370 
-880 



i7o: 



-815 
-630 



Table 5: Peculiar velocities for A2163 and A1689 and 
all considered contributions to the uncertainty. 



8. Conclusion 

The kinematic S-Z effect has the potential to be a 
singularly sensitive probe of the peculiar velocities of 
distant galaxy clusters. The realization of this poten- 
tial requires the separation of the thermal and kine- 
matic components of the S-Z effect. This can only 
be practically carried out by observations at mm and 
sub-mm wavelengths. We have used observations of 
the S-Z effect in three mm- wavelength bands to place 
limits on the peculiar velocities of two distant clus- 
ters. 

The cluster A2163 was observed in the 2.1, 1.4, 
and 1.1 mm bands. The 2.1 mm band observations 
are used to determine the position and morphology 
of the peak S-Z surface brightness. The positions 
of the X-ray and S-Z peak surface brightnesses are 
consistent when uncertainties in the ROSAT/PSPC 
and SuZIE pointing are included. The S-Z surface 
brightness is consistent with the X-ray derived den- 
sity model and the assumption of an isothermal IC 
gas. From the 1.1mm band observations, we report 
a significant detection of the S-Z thermal effect as an 



increment. Simultaneously fitting the 2.1, 1.4, and 
1.1 mm data with the isothermal kinematic and ther- 
mal component surface brightness models, we deter- 
mine yo = 3.62±.49x 10"-* and Vr = -H490;^3°kms"^ 
at 68.3% confidence. The uncertainties in Vr due to 
the thermal structure of the IC gas are found to be 
much smaller than the statistical uncertainty. 

The cluster A1689 was observed in the 2.1 and 
1.4 mm bands. The 2.1 mm band observations were 
used to determine the position of the peak S-Z surface 
brightness. The PSPC and S-Z peak surface bright- 
nesses are nearly coincident. Fits to the S-Z surface 
brightness indicate the observed morphology is con- 
sistent with the density model determined from the 
X-ray observations and isothermal IC gas. Simultane- 
ously fitting to the 2.1 and 1.4 mm data, we determine 
2/0 = 3.43 ± .59 X lO"'' and Vr = +170t5TOkms~^ at 
68% confidence. 

Unlike the determination of Hq, the determination 
of cluster peculiar velocities does not place high de- 
mands on the quality of the X-ray data. The results 
are constant over the allowed range of X-ray derived 
density models. For both clusters, we have deter- 
mined Vr over a broad range of X-ray temperature; 
the value of the peculiar velocity depends weakly on 
the assumed Tg. Astrophysical confusion of the pe- 
culiar velocity measurements in the two clusters is 
likely to be dominated by primary anisotropies of 
the CMB at a level < 300 kms"^ All the sources 
of uncertainty to the determination of cluster pecu- 
liar velocities, that we have considered, are listed in 
Table pi Including uncertainties due to the uncer- 
tainty in the IC gas temperature, thermal structure, 
density model, calibration, and astrophysical confu- 
-t-490J:88o''kms"^ for A2163, and 
-^ for A1689, both at 68% confi- 
dence. 

Even though we have not yet reached sufficient sen- 
sitivity to determine individual peculiar velocities, it 
is possible to use results like those presented here, to 
probe the velocity field on very large scales. For exam- 
ple, Lauer & Postman (1994) (LP) report a bulk flow 
of galaxy clusters within a sphere of 15, 000 kms~ of 
VLp = -1-730 ± 174 kms"^ A2163 and A1689 lie along 
lines of sight separated from the direction of the LP 
bulk flow by only w 29° and Ri 17° respectively. If 
we assume that all the radial peculiar velocity of the 
clusters is due to a bulk flow in the direction of the 
LP results, we find vtuik ~ +280tlfo'kms^^ at 68.3% 
confidence. The results of this work cannot rule out 



sion, we find v> 
Vr = -hl70 



630 ^'^^^° 
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the continuation of the LP bulk flow to z ^ .2. A 
modest increase in sensitivity coupled with the obser- 
vation of clusters in the opposite direction to the LP 
flow (as a systematic check) would make a definitive 
test possible. 

The sensitivity of the S-Z observations described 
here are limited by differential atmospheric emission. 
It should be possible to use the distinct spectral sig- 
natures of the S-Z effects to separate them from at- 
mospheric noise. This requires simultaneous measure- 
ments along a given line of sight in at least three mm- 
wavelength bands. We have constructed and begun 
observations with an instrument capable of such mea- 
surements. This instrument should eventually have 
the capability of determining cluster peculiar veloci- 
ties limited by the astrophysical confusion of the S-Z 
effect. 
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